Optical lattice quantum Hall effect 



Rebecca N. Palmer 1 ' 2 , Alexander Klein 1 ' 3 , and Dieter Jaksch 1,3 
1 Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom 
2 School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom 
3 Keble College, Parks Road, Oxford 0X1 SPG, United Kingdom 
(Dated: July 8, 2008) 

We explore the behavior of interacting bosonic atoms in an optical lattice subject to a large arti- 
ficial magnetic field. We extend earlier investigations of this system where the number of magnetic 
flux quanta per unit cell a is close to a simple rational number [Phys. Rev. Lett. 96, 180407 (2006)]. 
Interesting topological states such as the Laughlin and Read-Rezayi states can occur even if the 
atoms experience a weak trapping potential in one direction. An explicit numerical calculation 
near a = 1/2 shows that the system exhibits a striped vortex lattice phase of one species, which 
is analogous to the behavior of a iuio-species system for small a. We also investigate methods to 
probe the encountered states. These include spatial correlation functions and the measurement of 
noise correlations in time-of-flight expanded atomic clouds. Characteristic differences arise which 
allow for an identification of the respective quantum Hall states. We furthermore discuss that a 
counterintuitive flow of the Hall current occurs for certain values of a. 

PACS numbers: 37.10. Jk,03. 75. Lm,73.43.-f 



I. INTRODUCTION 

Since their discovery, the integer [l[ and fractional 
quantum Hall effects [3, Q have attracted much atten- 
tion on the theoretical as well as on the experimental 
side. They were first observed in semi-conductor samples 
which confine the electrons in a quasi two-dimensional 
layer and are characterized by the ratio v of particles 
to magnetic flux quanta penetrating the layer. For the 
integer quantum Hall effect, this ratio is an integer num- 
ber, whereas the fractional quantum Hall effect exhibits 
a simple rational (non-integer) v and has interparticle in- 
teractions as an essential component. Some of the excita- 
tions of the fractional quantum Hall effect are promising 
candidates for topological quantum computing [4|, since 
they are separated from the ground state by an energy 
gap and they might exhibit fractional, anyonic statistics 

la a. 

Despite the work which has been put into this field, the 
mechanisms which lead to the Hall effects are not com- 
pletely understood 0, HI • Especially theoretical investi- 
gations are hampered by the complexity of the system: 
A many-body quantum description is necessary to cap- 
ture the relevant effects such as the strong correlation of 
the electrons, which means that the calculational effort 
scales exponentially with the number of particles to be 
described. To gain insight into the behavior of the quan- 
tum Hall effects it is therefore worthwhile to investigate 
whether these phenomena occur in alternative systems 
as well. 

Ultracold atoms in optical lattices are a system where 
the quantum Hall effect can indeed be observed. As some 
of the present authors have shown in an earlier publica- 
tion |9j, for a suitable choice of the external parameters 
various topological states occur in these systems such as 
the Laughlin state or the Read Rezayi states [TTI | . 
These states have been derived by developing a model 



that is valid near simple rational numbers of magnetic 
flux quanta per lattice cell, a. We showed that the oc- 
curring states can be distinguished by shot noise and 
Hall current measurements, where in the latter case un- 
expected sign changes can occur. 

In the present paper we extend these investigations. 
We will give details of how to expand the Hamiltonian 
describing cold bosonic atoms in an optical lattice subject 
to an artificial magnetic field for the case that a is close 
to a simple rational a c = l/n with I and n small integers. 
We will also give more insight into how to distinguish 
the occurring states by noise correlation or Hall current 
measurements. In addition to these extensions, we focus 
on other phenomena occurring in our setup. By using 
a Gutzwiller ansatz we demonstrate that for a close to 
a c = 1/2 and large filling fractions a two-component-like 
striped vortex state occurs which shows that the analyti- 
cal description derived in Q is meaningful. The presence 
of the magnetic field leads to an additional small scale 
structure for values of a close to a c , which can be made 
visible in time-of-flight expansions as we will illustrate 
fora c = l/2. 

The applicability of our scheme is ensured by the ex- 
perimental and theoretical successes in the field of optical 
lattices during the last couple of years. Theoretical in- 
vestigations have shown that a whole wealth of models 
can be simulated [2] , such as the Hubbard Hamiltonian 
for bosons, fermions, or mixtures (l3l. [Til], spin-spin in- 
teractions [HI, [l(| , high-temperature sup erconductivity 
(Hill, or the formation of polarons [H, H3, EH ■ Theo- 
retical investigations show that these models exhibit rich 
phase diagrams which should be accessible with current 
experimental techniques [TH, [22], H3, [24], The nec- 

essary low temperatures can be achieved by appropriate 
cooling methods [H, [U H, H^j. It is also possible to 
implement artificial magnetic fields in this setup, for ex- 
ample by rotating the lattice [3(| HH, H3, HU . This comes, 
however, with the disadvantage that a centrifugal term 
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occurs, which has to be balanced by an external trap- 
ping potential. To circumvent this problem alternative 
schemes have been developed, exploiting Raman assisted 
hopping [34], HH, oscillating quadrupole fields [36^ 1 , o r 
laser beams with orbital angular momentum [37], [38|, l39( . 
Some of these schemes even exhibit non-Abelian gauge 
potentials [Z(J, Sl[ , allowing for the investigation of new 
phenomena in the integer quantum Hall effect [42| or 
the observation of the non-Abelian Aharanov-Bohm ef- 
fect m. 

In the presence of a lattice, the strength of the artificial 
magnetic field is suitably characterized by the number 
of effective magnetic flux quanta per unit cell, a, which 
for a conventional semiconductor setup is typically very 
small, a < 1, In the optical lattice setup, however, it 
is possible to achieve values close to a = 1/2, a regime 
which is almost impossible to achieve in the condensed 
matter case due to the required huge magnetic fields. 
This allows for the experimental investigation of so far 
unobserved phenomena such as the Hofstadter butterfly 
[34L [44|, |45|, |4(| , a fractal energy spectrum originally pre- 
dicted for electrons subject to a large magnetic field. It 
is also worthwhile to investigate the energy gap of any- 
onic excitations near a — 1/2. If it is sufficiently large, it 
would make topological quantum computing more feasi- 
ble since thermally excited anyons are the major source 
of errors in these schemes (3| ■ 

There have been previous investigations of ultracold 
atoms in optical lattices subject to large artificial mag- 
netic fields or rotation. For instance, in [53] a rotation 
scheme of an optical lattice was proposed which leads 
to quantum Hall and spin liquid properties of a spin-1 
boson cluster. The dependence of the state on the rota- 
tion frequency has been investigated in Ref. [4q . where 
it was shown that one can change the state adiabati- 
cally from a Gaussian to a Laughlin state. Transport 
properties of one and two particles in a rotating lattice 
have been investigated in Ref. (4|| using a Kubo formal- 
ism. Recently, the influence of boundary conditions on 
the Hall states was investigated and Chern numbers were 
employed to characterize the signature of the topologi- 
cal states [50|, HH ■ These authors also discussed that the 
quantum Hall states can be made more stable and the en- 
ergy gap to excitations can be increased by using atoms 
with a strong dipole-dipole interaction such as chromium 

An experimental implementation of the proposed 
methods and the observation of the predicted effects 
should be possible in the near future [lj|. The creation 
of fast rotating traps has already been demonstrated in 
Refs. [H, HII' which have allowed one to investigate ef- 
fects such as vortex pinning [54| . Loading an optical 
lattice with a well-defined number of atoms, thereby al- 
lowing one to tune v independently from the strength 
of the artificial magnetic field, has been achieved by ex- 
ploiting spin oscillations [H, [B(| [57]], and the number 
statistics of ultracold atom systems has been examined 
in [58l . [59! ] . Also, measurement schemes have been im- 



plemented which allow for an investigation of second or- 
der correlation functions such as noise correlations [6C| . 
thereby revealing additional quantum properties of the 
states. 

Our paper is organized as follows. In Sec. HH we in- 
troduce the Hamiltonian which describes the atoms in 
the optical lattice subject to an artificial magnetic field 
in Landau gauge. For weak magnetic fields, a continuum 
approximation is applied as shown in Sec. lIIIl We extend 
our investigations in Sec. lIVI to higher artificial magnetic 
fields, such that a is close to simple rational values where 
the system is well-described by multi-component contin- 
uous wave functions. In Sec. [V] we present how the differ- 
ent quantum Hall states can be distinguished from each 
other. The used methods include time-of-flight expan- 
sions, measuring noise correlations in the time-of-flight 
expanded clouds, and observing the occurring Hall mass 
current. We conclude in Sec. IVII 



II. HAMILTONIAN 

Throughout this work we assume a two-dimensional 
optical lattice [H 03, which can be experimentally 
achieved for example by using a three-dimensional lattice 
with the hopping into one direction strongly suppressed 
[6l| . The atoms are then confined in two-dimensional 
planes, of which we consider only one in the follow- 
ing. The artificial magnetic field term is created using 
one of the several methods discussed in the Introduction 
[34l [35l [3^ . l37l [H, HH. We furthermore assume that 
the total potential energy experienced by the atoms is 
given by V(p, g), where p, q label the lattice sites in the 
x- and y-directions, respectively. In Landau gauge, which 
is more convenient for our purposes, the Hamiltonian de- 
scribing the atoms is given by [H, O, [36| 

H = - ^2 (J x e 2maq al q a p -i. q + J y al q a p , q -i + H.c.) 
+ V(p, q)al q a p , q + g-a^aJ.gOp^Op.g , 

(1) 

where the magnetic field strength is parameterized by the 
number of flux quanta per lattice cell a. The operator 
a p q creates an atom in lattice site (p, q), J x and J y are the 
hopping constants in the x and y directions, respectively, 
and U describes the strength of the on-site interaction 
between the atoms. Due to the artificial magnetic field 
terms, the atoms acquire a phase 2na when they hop 
around a lattice cell once. Since the value of a is only 
defined modulus one, we may restrict it to < a < 1. 

It is instructive to compare our setup to other sys- 
tems in which the quantum Hall effect can be observed, 
such as a semiconductor structure or a rotating ultracold 
atomic gas without an optical lattice. The main charac- 
teristics and typical values are shown in Table [T] Due to 
the small lattice spacing d ~ 10~ 10 m in the original solid 
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state semiconductor systems, for instance, Si or GaAs, 
typically only small values a — eBd 2 /2hn ~ 1CP 4 can 
be achieved with available magnetic fields B. Here, e is 
the charge of the electron. Although this can be over- 
come by using superlattice structures [H, [62j , defects or 
impurities are inevitable in a real crystal. In contrast, 
the optical lattice setup allows for a virtually defect-free 
implementation of the quantum Hall Hamiltonian, and 
values up to a ~ 1 can be achieved [H, [H, [H, [50| . 
Comparing the cyclotron frequency fl — iraH/fhd 2 (for 
a definition of the effective mass fh in the lattice setup 
see the next section) to the temperatures which can be 
achieved in the respective realizations shows that simi- 
lar regimes can be realized with all three methods when 
appropriate coolin g te chniques are exploited in the atom 
setup [26|, [27], [H, |29| . Low ratios of temperature over 
magnetic field are important to avoid thermal excita- 
tions which would spoil the applicability of the system 
for topological quantum computation [7J. 

In contrast to the semiconductor scheme, the two 
atomic realizations provide more opportunities for tun- 
ability and probing the states. Both atom setups al- 
low for additional measurements to be carried out, such 
as time-of-flight expansions and measuring second order 
correlation functions [6(J US] • The lattice setup further- 
more enables an easy tuning of the system parameters 
such as the hopping or the interaction between the atoms 
[3]. In a standard optical lattice setup, this interac- 
tion is given by a contact term S(r). By using certain 
atomic species such as chromium one can also implement 
dipole-dipole interactions which might stabilize the quan- 
tum Hall states [50j . It is also possible to extend the 
contact interaction by immersing the optical lattice into 
a Bose-Einstein condensate, where an additional off-site 
interaction occurs [13, ufl, k5J| . 



Compared to the rotating setup, in a lattice there ex- 
ist methods which do not involve a centrifugal term that 
has to be balanced by an additional external potential. 
Although experiments with rotating traps already exist 
[52I l53l . [Hij], reaching the quantum Hall regime is chal- 
lenging due to the need of balancing the centrifugal force. 



III. SMALL a LIMIT: CONTINUUM 
APPROXIMATION 



In earlier publications 0, [36[, it was shown that for 
small values of a the influence of the optical lattice is 
negligible and we can approximate the lattice gas by a 
continuous wave function. In the present section, we give 
a more detailed account of this analysis. We first consider 
a single particle in the optical lattice described by the 
discrete wave function | t/>) = ^2 p ip(p, q)a pq | 0), which 
obeys the normalization condition \ip(p, q)\ 2 = L 



The Hamiltonian Eq. (jTJ) acts on this wave function as 
H\i>) = J2{-U^ iag ^(p -!,«) + e- 2li °V(p + 1, <?)] 



J v [il>(p,q-l)+1>faq+l)] 
■V(p,q)1>(p,q)}al q \0). 



(2) 



For a <C 1 and a weak trapping potential V(p,q), 
the wave function varies only slowly from one lattice 
site to the neighboring ones, and we can approximate 
the state by a continuous wave function <p(x,y), where 
ip(p,q) — d<ft(pd,qd) with d the lattice spacing and 
J\(p(x, y)\ 2 dxdy = 1. The dynamics of this wave function 
is governed by the Hamiltonian 



Ho — — J x 



J 11 



d 2 ( d wfr V 
h 2 \ dx d 2 J 



d2 (■* 9 



(3) 



V(x,y) 



This Hamiltonian can be transformed into a more fa- 
miliar form by defining the effective masses m xy = 
H 2 /2 J x , y d 2 , fh — ^Jm x m y , and the cyclotron frequency 
fl = ■nah/fhd 2 . After discarding a constant energy term 
we get 

This is the familiar single particle quantum Hall Hamil- 
tonian with an artificial "electric" field potential V(x,y) 
and anisotropic mass, which can be redefined into an 
anisotropic magnetic length. For deriving the energy lev- 
els we assume that the potential V(x, y) is constant in the 
x direction and forms a harmonic potential in the y di- 
rection, V(x, y) — m y u> 2 y 2 Jl. The Hamiltonian is then 
translational invariant in the x direction, which justifies 
the ansatz cj>(x,y) = exTp(—iKx)F(y). The Hamiltonian 
Hq acts on this wave function as 



H n 6 = 



2 



2m y dy : 



+ \ m yU 2 cS {y ■ 



lo 2 K 2 K 2 



where to e s — vSH 2 



and y c — 2ilhK/u> 2 s f 



(5) 
This 



expression describes a displaced harmonic oscillator with 
an energy offset depending on K. For u) = we retain the 
usual Landau levels, whereas for ui ^ the degeneracy 
within one Landau level is lifted and the energies are 
given by 



E = n LL 



u 2 h 2 K 2 
2m x uj 2 B 



(6) 



where till > is an integer. 

The continuum approximation can be extended to the 
case of more than one particle in the lattice, provided 
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System 


Electrons (GaAs) 


Rotating gas (Rb) 


Lattice (Rb) 


Particles 


Electrons 


Atoms 


Atoms 


Statistics 


Fermi 


Bose 


Bose 


Interaction potential 


1/r 


8(r) 


S(r) 


Confinement 


Sharp edges 


Smooth trap 


Smooth trap 


Physical parameters 








System diameter 


~ 10" 2 m 


10 -5 - 10 _4 m 


10 -5 - 10 _4 m 


Number of particles N 


~10 n 


10 5 - 10 6 


~ 10 4 


2D number density q 


1 - 2 x 10 15 m~ 2 


~ 10 12 m" 2 


10 10 - 10 n m- 2 


Cyclotron frequency (CF) Q, 


~ 10 14 Hz 


10 2 - 10 3 Hz 


~ 10 3 Hz 


Free mass mo 


9 x 10~ 31 kg 


1.5 x 10~ 25 kg 


1.5 x 10" 25 kg 


Effective mass m 


6 x 10~ 32 kg 


1.5 x 10" 25 kg 


l -24 _ lQ-23 kg 


Lattice spacing d 


2 x 10 _1() m 


n/a 


~ 10" 6 m 


Temperature T 


0.1 -0.3K 


~ 10 _9 K 


~ 10~ U K 


Dimensionless parameters 








a (flux quanta/unit cell) 


ht 4 


n/a 


0.1-0.5 


v (particles/flux quanta) 


0.2 - 10 


10 3 - 10 4 


1 - 10 


Qu/hil (interaction/CF) 


n/a 


0.1 - 1 


~ 0.1 


k B T/hn (temperature/CF) 


10" 3 


10" 1 - 1 


~ 10~ 3 


Internal states 


spin or bilayer 


hyperfine 


hyperfine 


Tunability 








n, e ,T 


yes 


yes 


yes 


um 


no 


yes 


yes 


m x /m y 


no 


no 


yes 



TABLE I: Comparison of quantum Hall systems. The regime of interest is v ~ 1, and hQ. ^> gu S> ksT. The values for the 
BEC are from current experiments [s3. l53|. which are still slightly outside this regime. The low temperatures for the lattice 
system can be achieved by using, for instance, the methods in Refs. [26l. |27|. |2S|. \2ty . For the lattice, the data shown is for a 
single plane only. If the two-dimensional plane is realized using a three-dimensional lattice with the hopping in one direction 
strongly suppressed, the quantum Hall effect occurs independently in each plane and multiple planes can be used to increase 
the signal strength without increasing v. 



that the wave function still varies slowly over the distance 
d between two lattice sites. To achieve this, the average 
distance between the particles needs to be large compared 
to d, which restricts the particle density to g -C 1/d 2 . In 
this case, the continuum Hamiltonian is given by 

H = Hoixj,Vj) + t;Y1 S ( Xi ~ x j) 5 (y* ~ Vi) ' ( 7 ) 

3 i,3 

where u = Ud 2 . This Hamiltonian describes the quantum 
Hall effect with a contact interaction term, in contrast 
to the typically screened 1/r interaction of the Coulomb 
potential for electrons in solids. In analogy to the solid 
state fractional quantum Hall effect we define the filling 
factor v = hgn/rhfl = gd 2 /a. Since, in this section, 
we are interested in the limits a <C 1 and g -C 1/d 2 we 
will hence find the same states as in a continuum bosonic 
quantum Hall system, for example, the v = 1/2 Laughlin 
state, the v = 1. 3/2, . . . Read-Rezayi states, or a vortex 
lattice dElil. 

The strength of the interparticle interaction u plays a 
crucial role if the setup is to be used for topological quan- 
tum computing. For this, anyonic quasiholes are created 



and they are moved around each other using focussed 
lasers, thereby inducing qubit operations [3l[. However, 
for the Read-Rezayi (non-Abelian anyon) states to ap- 
pear the interaction must be weak compared to the spac- 
ing between the Landau levels given by 2hQ, to avoid ex- 
citations to higher levels. This reduces the quasiparticle- 
quasihole pair creation gap A g ~ wmf2/27rft, where we 
assume m = fh = m x = m y . As thermally created 
anyons moving around the computational anyons are a 
source of error, such computation would require low tem- 
peratures T <C umfl/27rhkB <C 2hQ/kB, which for typi- 
cal experimental setups are on the order of a few lOnK. 
These temperatures might be reached using the methods 
of USES 11E HI- Alternatively, it has been suggested 
that non-Abelian field quantum Hall states [H, Ell may 
offer non-Abelian excitations in the lowest density state, 
allowing a strong interaction to be used to increase the 
gap to ~ Ml, but this has yet to be confirmed. 
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IV. NEAR SIMPLE RATIONAL a: or equivalently 

MULTI-COMPONENT WAVEFUNCTIONS 



We will now relax the condition of small a and investi- 
gate the properties of the quantum gas for an a close to 
simple rational values a ~ l/n, where I and n are small 
integers. The wave function can then be approximated 
by a set of n smooth, slowly varying functions which 
correspond to n different components of the gas. These 
components should not be confused with, for instance, 
different internal states of the atoms. They rather corre- 
spond to different small scale structures of the gas. Based 
on previous work [9J] , we will give details on how to derive 
these expressions and consider some special cases. We es- 
pecially calculate the ground state wave function for a set 
of representative parameters numerically, which confirms 
that for a close to 1/2 the two-component description is 
meaningful. For simplicity, in the following we will con- 
sider the case of isotropic hopping, that is, J x — J y = J, 
leading to effective masses m x = m y — m = m. 



> 



C 



(4C(a 



/3 2 ) 



>27r(a-a c ). (10) 



Here (3 = md ui/hir is the dimensionless trap strength. 
The first inequality shows that the range over which the 
multi-component ansatz is valid gets narrower for larger 
denominators n of a c . We also see that the presence of 
the trap is important to ensure the validity of the sec- 
ond inequality. A stronger trap makes it easier to fulfill 
the second condition, however, a too strong trap will ul- 
timately lead to a contradiction with the first condition. 
Indeed, in Ref. [9[ it was shown that for small n the 
above wave functions compare very well with numerical 
exact calculations and overlaps of more than 99% can be 
achieved for a~ a c and appropriately strong traps. 



B. Interacting particles 



A. Single particle states 



Eq. ([8]) is extended to the many-particle case as 



As discussed in [9j, numerical calculations for simple 
a = l/n and weak external potentials V suggest that 
the single particle ground state functions exhibit an n- 
site periodic pattern superimposed on a smooth large- 
scale variation. This motivates the representation ip(np+ 

hm+j) = dJ2kXk(d{np + i),d(nq + j))v^\ where \k is 
a continuous, slowly varying function and an n x n 
matrix describing the small scale structure of the atomic 
gas. We find by expansion about a c = l/n (see Appendix 
fAl for details) that there are n degenerate matrices v( fe ) 
of the form vj,^ = e 2 ^ tpk ^ n v q+ k, where v is a fixed n- 
component vector for each l,n and the subscript q + k 
wraps around mod n. Furthermore, the \k obey the 
condition 



c h 2 d 2 X k 

2m dy 2 



C_ 

2m 



2mf2y — ih 



d_ 

Ox 



= E- 



id 2 



Xk , 



Xk + V(x, y)xk 
(8) 



where = h(a — a c )ir / (md 2 ) and C,Eq depend only on 
l,n. This formula reduces to Eq. Q for a c — 0/1, and 
agrees well with numerical calculations near other simple 
a c , especially a c — 1/2 Q. 

For this procedure to be consistent, the length scale l x 
over which Xk varies, which for a harmonic trap V(x, y) — 
\muj 2 y 2 is given by l x = {H 2 C/{ACm 2 n 2 + m 2 ^ 2 )) 1 / 4 , 
must be large compared to the "small scale" periodicity 
nd, but small enough that 2mQl x d/h <ti 1 (cf. Appendix 
E}, that is 



1 1 2mdQ 
nd l x h 



(9) 



H 



Ch 2 d 2 
2m dy 2 



-u G k u k 2 ,k 3 Mx{ 1 (x,y)x{ 2 (x,y) 

ki,k 2 ,k 3 ,k 4 

E 



x Xk 3 (x,y)xkA x iy) + 



in 



d 2 ' 



(11) 



where G klM MM = J2j v j+ki^j+k 2 ^j+k 3 ^j+k i /n if ki + 
ki = &3 + ki mod n, and otherwise. The conservation 
mod n is due to the x quasimomentum 2nk /nd carried by 
v( fc ). To give an example, for a c — 1/2 we find Gim = 
G2222 = 3/2 and G1212 = G2121 = G1122 = G2211 = 1/2. 
Particles with the same k interact more strongly because 
their are peaked on the same sites. 

For a c — 1/2, a change of basis from \i,2 to x± = 
Xi ± ix2 makes this effective Hamiltonian analogous to 
a bilayer fractional quantum Hall system (66j . with x± 
being the two "layers" . However, the interaction ratio is 
1:2 with the "interlayer" interaction being the stronger 
one, while it is equal or weaker in most other realizations 
of multicomponcnt fractional quantum Hall states. 



C. Some special states 

For the case that a c = 1/2 the system is equivalent to a 
two-component gas as argued in the previous section. It 
is well-established that for an interacting two-component 
gas the lowest Landau level state with highest density 
and zero interaction energy is the so-called 221 state [67| 
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defined by 

0221 (zi, z N / 2 , wi,.. .,w N/2 ) 



1>J 



x exp 



l>3 



hi 



Wj 



M 2 /4 



V4 



(12) 



where there are N/2 particles in one component with co- 
ordinates z = (x+iy)/rQ and N/2 in the other component 
with coordinates w = (x + iy)/ro, and tq — [h/2mri] 1 ^ 2 . 
For the lattice setup at a = 1/2, the two components 
are replaced by the two "layers" , and the effective filling 
factor of this state is defined with respect to f2 and given 
by v = hgn/mfl — 2/3. The state can be extended to a 
general a c — l/n with v = n/(n + 1). 

Since this state has exactly zero interaction energy, 
it would be the lowest step of the density profile in a 
slowly varying external potential. Adding more particles 
to the system leads to a trade-off between an increasing 
potential energy when the particle is added to the "edge" 
of the 221 state and an increased interaction energy if 
the particle is added at the center of the trap where the 
potential energy is lowest. This eventually leads to a 
stepped density profile. So far it was not possible to 
determine with certainty what the next higher step looks 
like, but it might be either a v — 1 Read-Rezayi [llj state 
or a v = 4/3 non-Abelian spin singlet (NASS) state [68j |. 
Depending on the energies of the states, which are not 
exactly known, also both steps can occur. Higher states 
might be the Read-Rezayi states with v — 3/2, v = 2, 
and so on, eventually leading to a vortex lattice [6^| . 

It is still unknown where exactly the transition to this 
vortex lattice phase occurs. The Lindemann criterion 
[6^ . l70l ] gives the order of magnitude estimate v ~ 10 
for both one- and two-component systems, whereas nu- 
merical studies of the single component system [65j | give 
v ~ 2 — 6. We hence conjecture that our two- "layer" 
system has a similar critical v. To investigate the vor- 
tex lattice at a large filling factor and a close to 1/2, we 
have done numerical simulations using imaginary time 
propagation of the Gutzwiller ansatz 



iv>)=n 



1,3 



n=0 



(13) 



where the state is specified by the complex numbers Ci t j n . 
This ansatz can describe both superfluid vortex lattice 
and Mott insulator states, but in general not off-site cor- 
related states such as the fractional quantum Hall states. 

A typical density distribution of our calculations using 
a random initial state is shown in Fig. [TJa) . It exhibits 
alternating stripes as have been predicted for a contin- 
uum two-component system [72|, [73|. However, in our 
case the two components are not two different species or 
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FIG. 1: (Color online) Numerical approximation to the 
ground state of the ot c — 1/2 vortex lattice phase, calculated 
using imaginary time propagation of a Gutzwiller ansatz. (a) 
Density of the atom distribution, black stands for no den- 
sity, white means high superfluid density, red (hardly visible) 
means high Mott insulator density, (b) Population of the two 
"layers", black and white depict the two components, color 
encodes the phase of the superposition of those two compo- 
nents. The parameters are a = 0.48, U = 0.1J, and v = 21.2. 



states of one species, but represent the two "layers" cor- 
responding to the two components of the wave function. 
For a better distinction these two components are plotted 
with different colors in Fig. [TJb) . A direct measurement 
of these two components analogous to the measurement 
of two different species is not possible. Instead, to make 
the two components visible we calculated the overlaps of 
the small scale solutions v( fe ) with groups of lattice sites 
of the Gutzwiller solution. These overlaps correspond to 
the large-scale functions Xk at the respective lattice sites. 
Our calculations show that the densities described by the 
two components indeed add up to the total density of the 
system and confirm that at least for the chosen param- 
eters the two-component description is meaningful. We 
find that the stripes straighten out only very slowly with 
increasing imaginary time (or decreasing temperature), 
similar to the findings in [73j, with a slight tendency to 
phase separation. The perfectly straight stripes observed 
in [jl] exist only in the zero temperature limit. At ex- 
actly a — a c = l/n, we have v = oo and one has an n 
component superfluid with no vortices. This may be the 
possible superfluid phases found at a = 1/3,1/2 in the 
exact diagonalizations of S0rensen and coworkers [36j . 



V. PROBING THE STATES 

After having demonstrated that a very wide range of 
different states can occur within our setup, the ques- 
tion arises how these states can be distinguished from 
each other. Fractional quantum Hall states for which 
Eq. (|11[) is valid are characterized by a fixed and 
hence the density for these states will obey the condi- 
tion g a | a — a c \. For "ordinary" fractional quantum 
Hall states with a <^ 1 the characteristic parameter is is, 
which consequently leads to a density that obeys g oc a 
as discussed in Sec. IIII1 In contrast, for Mott insulat- 
ing states, the density is fixed to integer filling and does 
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neither depend on a nor on a — a c . An in-trap density 
profile measurement could hence distinguish these states 
by comparing measurements at different artificial mag- 
netic fields (and thus different a), as well as measuring 
the ratio v (or z>), at which these incompressible states 
occur. 

However, a more clear-cut distinction between the oc- 
curring states is possible by using methods such as time- 
of-flight expansions, and measuring noise correlations or 
mass currents. In the following sections we show that 
they reveal additional information which helps to iden- 
tify the respective states. We should note that some of 
our results rely on Monte-Carlo simulations of small sys- 
tems or rough estimates of particle energies and thus they 
should be considered as qualitative estimates only. 



after the release 



A. Time-of-flight expansion 



The time-of-flight expansion is a standard measure- 
ment tool for ultracold atomic systems. All potentials 
(including the artificial magnetic field) are turned off in- 
stantaneously and the atomic cloud is allowed to expand 
freely for a certain time before measuring its density dis- 
tribution. Since for an optical lattice setup interactions 
during the expansion can be neglected and the time can 
be chosen long enough such that the cloud expands to 
several times its initial size the result of this process is 
[zl 



£ m (X) = \w(Xm /M)\ 



2 E 

xi,x 2 



Qi (xi;x 2 )e 



iX- (xi — X2)rno /hi 



Here w is the Fourier transform of the Wannier func- 
tion describing an atom in a single lattice site, Q\ is the 
one-particle density matrix, the sum is taken over all lat- 
tice site vectors Xj, g m (X.) is the measured density at 
position X and time t after release, and mo is the free 
mass. Note that the "momentum" measured by this pro- 
cess is the free space momentum —ih\7. It is not to be 
confused with the momentum of the system Hamiltonian 
{—ihd x — 2m51y, —ihd y ). Furthermore, the operators &j 
need to be transformed into the laboratory frame, and 
Qi is expressed in terms of these laboratory frame opera- 
tors. This is because the sudden switch-off of the artificial 
magnetic field does not obey Maxwell's equations, which 
breaks gauge invariance. 

As shown in [75], for a low a one can calculate the 
single particle density matrix q\ of the constant density 
lowest Landau level states, which include the Laughlin, 
Read-Rezayi, 221, and NASS states. For large a, similar 
expressions with an additional factor ~J2 k v( fc )*v( fc ) from 
the small scale structure apply. This gives for the density 



g m (X) = |*(Xm /M)| 2 ^^v'W*(Ov'W(j)-^ 



x exp 





i,j k 
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' Xmn 


y% + Vj 






'Ymo 




ht 


Arl 



2irnrl 



[Xi Xj ) 



( W - y s ) , (15) 



where the vectors Xj = {xj,yj) label the lattice sites, 
and the prime indicates that we have taken the v'^ in 
symmetric gauge, which doubles their period to 2n. Thus 
the arguments of the v'( fc ) have to be taken modulus 2n. 

Summing first over Xi — x 2 and approximating this 
by an integral, we see that each X! + x 2 contributes 
an approximate Gaussian centered on (— y\ — y2,xi + 
X2)ht/Amorl, and then summing over xi +x 2 we get the 
in-trap density profile smeared out on a scale ro , rotated 
through a right angle and scaled by ht/Am^rl. Physically 
this happens because lowest Landau level wavefunctions 
are very similar to rigid body rotation, and once released 
from the field the atoms fly apart as if they were rotating. 
Because the sum is in fact a discrete one, the pattern re- 
peats every reciprocal lattice cell (size 2-rtht / dm ) with a 
decaying envelope w, and high a states can have multi- 
ple peaks per reciprocal lattice cell from the small scale 
structure. 

Fig. [2] shows some numerically calculated examples, 
the pattern repeats every reciprocal lattice cell with a 
slow decay given by w [76|. As expected from the above, 
fractional quantum Hall states have stepped peaks, fur- 
ther distinguishing them from the featureless w of a Mott 
insulator [6fJ, while vortex lattice states have Thomas- 
Fermi (inverted parabola) peaks. Low a states have one 
peak per reciprocal lattice cell while awl/2 states have 
4. The extra peaks (previously described for the vortex 
lattice in [33j]) confirm the existence of a small scale 
structure. However, they do not show whether one or 
more "layers" are occupied for multi- "layer" systems. 

As a time-of-flight expansion only measures single par- 
ticle properties, it cannot detect correlations, so different 
fractional quantum Hall states are distinguished only by 
their density. In the next section, we will consider an 
enhancement of this method that can measure the corre- 
lations of the states as well. 



B. Correlation measurement 

The two-particle correlations of a quantum state man- 
ifest themselves in shot noise correlations of the de nsity 
distribution after a time-of-flight measurement (63l . l77l |. 
as was recently demonstrated experimentally for the 
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FIG. 2: Numerically calculated time-of- flight expansions in 
symmetric gauge of (a,c) a Laughlin and Read-Rezayi state 
= 1/2, 1), and (b,d) a 221 and NASS state (u = 2/3,4/3). 
The initial states are circular with diameter 100 lattice sites 
for the outer ring and 50 for the inner one. The artificial 
magnetic field is chosen such that (a) a = 0.005, (b) a = 
0.505, (c) a = 0.01, and (d) a = 0.51. The resulting larger 
fields Q or Q, in (c,d), respectively, lead to a wider expansion 
of the cloud. The shown expansions do not include the w slow 
decay, as the speed of this is implementation dependent. A 
darker color corresponds to a higher density of particles, and 
figures (a) and (b) are to the same brightness scale, as are (c) 
and (d). Shown is the first Brillouin zone, and K = mo'K/ht. 



Mott insulator [60( ■ At time t after release from an opti- 
cal lattice these correlations are given by [H[ 

G(r,r') cxi- 6 £ ^'■^+^n'^ r '\a\a)a r <x v ) 

ii'jji 

-(n(v))(h(r')), 

(16) 

where i, j, i' , f run over all lattice sites, R^/ is the dis- 
placement vector from site i to site i , and Q(r) = mor/fti 
with too the free mass. For states described by a contin- 
uum wavefunction the expectation values over the quan- 
tum operators are given by (ajajaj'fij/) = d 2 £> 2 (i, j; i', j')- 
Here, Q2(i,j',i',j') is the continuum two-particle density 
matrix, which for lowest Landau level states with con- 
stant density is of the form [75j 



e 2 (i,j;i',j')=e-* [l * i|2+l ' 



' + +1*4 



-2(z' Zi ,+z' Zj ,)] 



27rnr 2 



]((Zi - Zj)* (Zi> -Zj-)). 



(17) 

Note that v = v and n = 1 for a c = 0. The function 
g(z) for arbitrary complex z is derived by analytic con- 




FIG. 3: Two-point functions for continuum fractional quan- 
tum Hall states: v = 1/2 Laughlin state (solid), v — 2/3 221 
state with gn — 322 (dashed) and gis = 1721 (dotted), and 
v = 1 Read-Rezayi state (dash-dotted). 



tinuation of the ordinary two-point correlation function 



g(\zi 



If the behavior of a two-component state 



close to a = 1/2 is to be investigated, we need to replace 

{a\a\ai>aji) =d 2 ^ Q 2 kl ' k2 \zi, zf, z v , z?) 



x v (fel) *(i)v (fcl) (0v (fe2) *(i)v (fe2) (i') , 

(18) 

where Q 2 — i s the continuum density matrix for 



two particles of the same type and g 2 



(1,2) 



Q^" 1 ^ for two 

particles of opposite type. The density matrices g 2 kl ' k ^ 
are described by expressions analogous to Eq. (| 1 T[) , where 
the function g(|z| 2 ) is replaced by gk ± ,k 2 (\ z \ 2 )- To obtain 
9kiM we St the series 



2ii(M 2 



922{\Z\ ) 
_ -W 2 /2 . 



E 

m>0 even 



ffi 2 (M 2 ) = 52i(M 2 ) = i + 



E 

m>0 



2c m 
to! 

2 Cm 

m! 



-M 2 /4 

1 

(19) 

-| z | 2 /4 

(20) 



to Monte Carlo data for the respective state. This works 
for any lowest Landau level state as shown in [jl, [79l |. 
where we take a sign change into account because our 
particles are bosons. Our calculations show that the 
Laughlin and 221 states have g(Q) = (see Fig. [3]), mean- 
ing that particles cannot come together and the interac- 
tion energy is hence zero, while the NASS and Read- 
Rezayi states have g(0) > and hence non-zero inter- 
action energy (the NASS state is not shown since it is 
qualitatively similar to the Read-Rezayi state). 

Eq. (|16p involves a sum over four complex variables, 
which makes its direct numerical evaluation computa- 



tionally intensive. To simplify this problem, we inte- 
grate over r + r', which gives a delta function setting 
Zi + Zj — Zi' — Zj> = 0, and removes all explicit depen- 
dence on z.i + Zj + Zi> + Zji . The resulting correlation, 
which (ignoring the t~ 6 prefactor) is a function of mo- 
mentum difference Ak = Q(r) — Q(r'), can hence be 
evaluated for an infinitely extended state by summing 
over just two complex variables, Zj — Zj and zy — Zji . 
However, this does not work for a finite size state be- 




cause the ranges of Zj — Zj and zy — Zj> then do depend 
on Zi + Zj + Zy + Zji , and assuming an infinite size state 
while necessarily summing over a finite range introduces 
the possibility of numerical artifacts. 

Some numerically calculated examples of the noise cor- 
relations are shown in Fig. [4J where again one recipro- 
cal lattice cell is shown and the pattern repeats. Fully 
Bose condensed states, including ideal vortex lattices, 
have zero correlation. The Laughlin and 221 states are 
found to show near 100% anticorrelation at small Afc. 
Here we should stress that as this measurement works 
in Fourier space, this anticorrelation is not the case for 
all fully real-space anticorrelated states: the Mott insula- 
tor, for example, shows positive correlation (60l |. Higher 
density fractional quantum Hall states also show anti- 
correlation but of reduced strength. Our data suggests 
that Read-Rezayi states have a ringed pattern with their 
strongest anticorrelation (ss 50% for v = 1 and « 35% 
for v = 3/2) at a nonzero Afc, while the v = 4/3 NASS 
state has « 40% anticorrelation at zero Ak and no ring. 
The presence or absence of this ring might be used to 
distinguish between the Read-Rezayi and NASS states. 
Apart from this we observe that the small-scale struc- 
ture for a w 1/2 is also visible in the noise-correlations, 
see Figs. |4|d)-|4|f): Additional anticorrelation dips at 
Afc = (±7r/2d,±7r/2d), (±7r/2d,0), (0,±7r/2d) occur, 
which have a similar structure to the central ones at 
Afc = (0,0). 



"° 
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FIG. 4: Numerically calculated shot noise correlations for an 
infinite system in symmetric gauge. The states are (a) Laugh- 
lin, (b) v = 1 Read-Rezayi, (c) v = 3/2 Read-Rezayi, (d) 221, 
(e) v = 4/3 NASS, and (f) v = 3/2 Read-Rezayi. The arti- 
ficial magnetic field was chosen to give a — 0.01 in (a)-(c), 
and a — 0.51 in (d)-(f). Gray corresponds to no correlation, 
whereas black corresponds to 100% anticorrelation. We can- 
not rule out that the four dips at Afc = (±7r/4d, ±7r/4d) are 
artifacts arising from our approximations. 



C. Hall current and disorder 

So far, we discussed only the static properties of the 
fractional quantum Hall states and how to detect them. 
However, in solid state physics one of the most important 
observables is a dynamical property, namely the Hall cur- 
rent. In this section, we will therefore study the atomic 
case subject to a linear tilt and to disorder. 

When a linear potential V(x,y) = may is applied to 
a continuum fractional quantum Hall system, all states 
acquire a velocity a/ (20) at right angles to the potential 
gradient causing the Hall current. In a perfect continuum 
system this is an exact result, unaffected by interaction, 
because such an acceleration term is exactly canceled out 
by Galilean transforming Eq. |7]) to a reference frame 
moving with this velocity. In a lattice system the result 
is valid at sufficiently low a for the continuum approx- 
imation to apply; at higher a the lattice, which defines 
a rest frame, becomes important and a different veloc- 



ity can occur [9(. In particular, near simple rational a 
where Eq. (fTTj) is valid, the velocity is a/(20), which can 
be very different: for a < a c it even has the opposite 
sign, meaning that the current flows in the opposite di- 
rection. A negative Hall current was also observed for 
a single particle in a lattice using Kubo response theory 

For nonlinear potentials of large length scale compared 
to the magnetic length ro, the single particle eigenstates 
lie along the equipotential lines of the potential and the 
Hall current flows along those lines at the velocity given 
by the local potential gradient. This motion is not visible 
in equilibrium as the steps in the density profile lie along 
equipotentials as well. However, it can be made visible 
by putting the system out of equilibrium, for example, by 
suddenly changing the trapping potential. Relaxation to 
equilibrium will be slow because nonforward scattering 
is energetically forbidden in the fractional quantum Hall 
system, making the Hall current a supercurrent [80 . |81| . 

It has been shown that when sufficiently mild disorder 
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is added to a fractional quantum Hall system, some of 
the particles become localized and cannot carry current, 
but those which remain free move faster and the average 
velocity is still a/(2J7) [80L fsTI] . For example, for smooth 
disorder V(x, y) = may + V r (x, y), where V r is a random 
potential with length scale much larger than tq and zero 
average, the expression for the x direction of the velocity, 
v x = (dV/dy)/(2mCl), remains valid. Hence particles on 
closed equipotentials are confined to those lines, but the 
random part of v x averages to zero so the average velocity 
in the x direction remains a/ (20,). The current is hence 
determined by the width over which each extended level 
is occupied. A simple model of this is to describe the 
disorder by a density of states Qd(E), given by the num- 
ber of lowest Landau level states per unit energy interval 
and per unit area. This gives J Qd(E)dE — mCl/(%K). 
In the case of smooth disorder, Qd(E) is proportional to 
the probability distribution of the noise potential. For a 
linear geometry, we then have 

N = L fdy f dE - "j-i)Qd{E - f*j) > (21) 

J J j>0 

(22) 

where L is the length of the system along the x axis, I 
is the net current in the x direction, N the number of 
particles, and \i the chemical potential. The fractional 
quantum Hall states are taken at filling factors Uj and 
chemical potential uj, and V is suddenly changed from 
V\(y) to V%(y) to achieve the nonequilibrium situation. 
The exact shape of the disorder and thus of V\(y) and 
Vi(y) are not important for our qualitative investigations. 
These equations are valid for weak trapping potentials, 
where dV/dy is much smaller than the disorder term. 
Stronger potentials can break weakly localized states free, 
replacing the S(E — fij) density of extended states in I 
by a finite width distribution. 

In a square-well potential (as in a solid state system) 
the localized states create a finite range of filling factors 
over which a given number of fractional quantum Hall 
extended state levels are completely full, giving rise to the 
almost perfectly flat fractional quantum Hall plateaus. In 
a harmonic trap there will not be fractional quantum Hall 
plateaus, only corners each time a new extended level 
begins to fill, see Fig. \5\ Between those corners, there 
exist different fractional quantum Hall states. Unlike the 
square-well case, it is possible to obtain the complete 
distribution Qd(E) by measuring I against N (or fi). 

VI. CONCLUSIONS 

In our work we have studied the low and high field 
fractional quantum Hall effect in an optical lattice. For 
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FIG. 5: Dimensionless Hall current Id = 

I j (2am 2 \J 2uCl j h? j \S) against dimensionless number of 

atoms per unit length Nd = N /(2Lm\J 2ufl 3 / h 3 / lj) , where 
L is the length of the system, which is contained in a ID 
harmonic trap V(x,y) = muPy 2 /2 + yma at a rss 1/2. The 
linear term yma is turned on after the atoms have come to 
equilibrium in the trap, and the curves are for no disorder 
(dotted straight line), maximal disorder (constant density 
of states, black curve), and Lorentzian disorder of width 
(l/5)umQ/2irh (gray curve). The corresponding plot for 
small q is qualitatively similar Q]. 



small a, corresponding to low field strengths, the con- 
tinuum approximation can be employed to find the solu- 
tions of the Hamiltonian. For high fields corresponding 
to values of a close to simple rationals, a rs a c = l/n, 
we showed that the states can be approximated by an 
n-component wave function, where each component is 
described by a slowly varying, continuous function x ■ 
These functions are solutions of a differential equation 
similar to the a«l case, but with an artificial magnetic 
field term Q = fm(a — a c )/md 2 and an effective filling 
factor v = hgir/mfl instead of the usual fl and v. 

We have shown that many interesting topological 
states may occur in our setup, even if it is subject to a 
weak harmonic trapping potential in one direction. These 
states include the Laughlin, the Read-Rezayi, the 221, 
and the non-Abelian spin singlet (NASS) state. Numeri- 
cal investigations also showed the occurrence of a striped 
vortex lattice phase close to a w 1/2, as predicted for 
a two-component system at low a. In our case the two 
components did not correspond to two different atomic 
species, but to the two components of the wave func- 
tion. 

We have furthermore demonstrated that the above 
states can be distinguished by suitable measurements. 
Simple time-of-flight expansions reveal additional struc- 
tures for higher a states, but do not necessarily allow 
for distinguishing between the different quantum Hall 
states. More insight can be gained by measuring two- 
point correlation functions. For the Laughlin state, this 
function is zero for \z± — z>i\ =0, meaning that there are 
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never two atoms at the same place, whereas for the Read- 
Rezayi state with v = 1 the correlation function exhibits 
a finite value at zero distance. We, moreover, showed 
that the different states also exhibit different noise cor- 
relations, which can be employed to distinguish between 
them. Again, higher a states show additional structures 
due to their additional small scale symmetries. Further, 
our results suggest that the Read-Rezayi states exhibit a 
ring structure, which allows one to distinguish them from 
the Laughlin or 221 states. 

In analogy to the semiconductor Hall effects, it is possi- 
ble to displace the optical lattice and measure mass trans- 
port effects, which behave similarly to the condensed 
matter analogues. Especially Hall currents proportional 
to a — a c rather than a have been predicted, which may 
lead to a negative sign in the direction of the mass flow. 
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APPENDIX A: DERIVATION OF HIGH a 
STATES 

In this appendix we give details on the derivation of 
the Hamiltonian equation |(5J| for large a close to a ra- 
tional a c — l/n with small integers / and n. We assume 
that apart from the optical lattice the atomic gas ex- 
periences no trapping potential in the x direction, but 
might be trapped by a (slowly varying) harmonic poten- 
tial V(x, y) = \mijj 2 y 2 in the y-direction. Note that the 
mass in this potential is the effective mass m = h/2Jd 2 
introduced earlier. Motivated by previous findings [9(, we 
make for the wave function the ansatz ij){ns+i, nr + j) = 
d(j)j((nr + j)d)e lK ( ns+l ) d . Here, i, j = 1, ...,n and r, s are 
integers, and we define x — (ns + i)d, y — (nr + j)d. 
Since H is independent of x if the potential V is, the x 
quasimomentum K is conserved exactly for our choice of 
V . Let K — 2nkl/nd + K , where k is an integer and 
K <C l/n to yield an n site periodicity plus slow variation 
in the x direction. The action of Hamiltonian equation 
(HJ with U = gives 



H\4,) = -jdj2J2 



r,s i,j L 

+ h-linr + J - iyK{ns +l )d + i+1 ( nr + j + 1)e i^(n S +i)d 



Kr+j,ns+i I vaC ) 



(Al) 



V ( nr + US + ^nr+j.ns+i I vaC ) 



r,s i,j 



Using the fact that <j>j is only slowly varying for constant index j [9( , we can approximate derivatives by appropriate 
discrete differences and make use of the fact 



t>j±i{y) | l d(j) j±1 {y) l d 2 cj) j± i(y) 1 ,.,,,„ 
± d—dy— + 2 dy 2 " ^M^r + 3± l)d] 



(A2) 



Collecting all terms with the same creation operators in Eq. (|A1[) we can derive a Hamiltonian for <f>j , acting as 



h 2 ( 2njl 2mflyd 

■ cos — ■ — I Kd 



id 2 



1 



Pj-i 



h 
id 2 c 



i n 2 

J 2 y V] 2m 



tf> j+1 , 1 86 



3+1 



ld 2 c 



d 2 d dy 2 dy 2 



d 2 d dy 2 dy 2 



(A3) 



+ 0(d) . 
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If the condition 2m(lyd/h <C 1 holds over the range where the wave function is appreciable, namely y ~ l x , 
apply the Taylor expansion to the cosine, which after collecting terms in orders of d yields 



H 



19 j 



2md 2 
2m 



K 



2 cos 



{2ir(j-k)l 



2mfly 
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h 2 
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(iv — )2sm 



Mj ~ k)l 



dy 



Mi - k)i 



ill 



id 2 c 



y+i 



2 dy 2 



15V 
2 dy 2 



0(d). 



dy 



(A4) 



Note that for n = 1 or 2 the odd powers in ci cancel by symmetry and we get an expansion in d 2 , but that for larger 
n the expansion is in d. Define $ — (</>i . . . 4> n ) T and expand in powers of d: $ = + d<&^ + d 2 ^ 2 ^ + 0(d 3 ). 
Similarly expand the energy E — Eo/d 2 + E\jd + E2 + 0(d). We define the matrix 



A 



/ 2cos2vr(l - k)l/n 1 

1 2cos2?r(2 - k)l/n 1 

1 







\ 



(A5) 



1 2cos2?r(n- k)l/n ) 



The 0(1/ d 2 ) terms then become A $ (0) = 
— 2mE &°) / h 2 , so $^ ^ is an eigenvector of A , 
with the ground state having the largest eigenvalue. 
Assuming nondegenerate eigenvalues, $(°)(y) is hence 
proportional to the same (normalized) eigenvector v(% 
and can depend on y only in overall magnitude, i.e., 

$ (0) (y) = x(2/)v (fc) . 

Since changing k is equivalent to changing the origin 
of j, the eigenvalues are the same for all k, with ground 
state eigenvectors — Wj-fc, where v is the normalized 
ground state eigenvector for k = and the subscript j — k 
wraps around modulus n. This gives n degenerate ground 
states k = 0, ...,n — 1, which are orthogonal because of 
their different K values. For simplicity we take k — in 
the remainder of this derivation. 



For the 0(l/d) terms, define 

/ 1 
-1 1 
0-10 



-1\ 



(A6) 



\ 1 ■•• -1 / 



A2 = 2diag[sin(27r//n), sin(47rZ/n), sin(2n7r//n)] , 

where diag denotes a diagonal matrix with the argument 
as the entries on the diagonal. Define furthermore wi = 



A]V, W2 = A2V. This yields 
(A + 2mEo/h 2 )^(y) = ^x(y> + gwj 



K- 




The left hand side of Eq. (|A8|) is orthogonal to v 
because Ao + 2mEo/h 2 annihilates v and is Hermi- 
tian, while Wi is orthogonal to v because Ai is an- 
tisymmetric and hence so is v T AiV, but the latter 
is a number so it can only be antisymmetric if it is 
zero. Hence a solution can only exist if W2 is also or- 
thogonal to v and E\ = 0; this is the case for all 
a c = 1/8,1/7,1/6,1/5,1/4,2/7,1/3,3/8,2/5,3/7,1/2 
within numerical accuracy Q , but we have not been able 
to prove that it is always the case. 
For the 0(1) terms define 



/ 1 

1 1 

1 

Vl ••• 



1\ 







'•• 1 
1 0/ 



(A9) 



A4 = diag[cos(27rZ/n), cos(47rZ/n), cos(2nirl / n)] 



(A10) 



13 



giving 
(2mE 2 



m u) y 



x(y)v 



( ZrnQy V 1 d 2 x • 



A 



2m-B 
"ft 5 " 



8jj"i + [ K -—) ^ 2 

A + ^)*( 2 )=0. 

(All) 

We note that for nondegenerate eigenvalues of A the ex- 
pressions (A + 2mE / h 2 )~ 1 Wj , j = 1,2, are well-defined 



since both wi and w 2 are orthogonal to v. Taking the 
scalar product of Eq. (|A11[) with v gives a harmonic os- 
cillator equation for x(y), 



E 2 X = - 



h 2 Ci d 2 x h 2 C 2 ( 2mVLy 



-K) x+V(x,y)x, 



2m dy 2 2m \ h 

(A12) 

where the dimensionless constants C\ = v T A3v/2 — 
wf (A + 2mE /ti 2 )~ 1 w 1 and C 2 = v T A 4 v - wj (A + 
2mEo/h 2 )~ 1 W2 depend only on I and n, and by chang- 
ing the roles of x and y by using a gauge transformation 
one can show that C\ = C 2 = C. This oscillator has 
mass m/Ci, frequency w e gf = (4CiC2f2 2 + Ci^ 2 ) 1 / 2 , and 
center y c = 2hC 2 Kn/(AC 2 mn 2 + muj 2 ). 
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